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Abstract 

In this work we introduce a completely general Chapman Enskog procedure in which we 
divide the local distribution into an isotropic distribution with anisotropic corrections. We 
obtain a recursion relation on all integrals of the distribution function required in the derivation 
of the moment equations. We obtain the hydrodynamic equations in terms only of the first few 
moments of the isotropic part of an arbitrary local distribution function. 

The incompressible limit of the equations is completely independent of the form of the 
isotropic part of the distribution, whereas the energy equation in the compressible case contains 
an additional contribution to the heat flux. This additional term was also found by Boghosian 
and by Potiguar and Costa [21 El m the derivation of the Navier Stokes equations for Tsallis 
thcrmostatistics, and is the only additional term allowed by the Curie principle. 

1 Introduction 

Boltzmann's two crowning achievements were his definition of the micro canonical entropy and 
Maxwell Boltzmann distribution, and his kinetic equation and mechanical basis for irreversibility 
of macroscopic non-equilibrium systems. Boltzmann's famous H-theorem gave the first connection 
of the second law of thermodynamics with underlying reversible mechanics. In this paper we shall 
revisit both of these contributions. Firstly we shall review the resolution of the apparent conflict 
between reversible microdynamics and the second law. We then turn to the main subject of this 
paper, which is the consideration of solutions to the Boltzmann equation when the local velocity 
distribution is not a Maxwell-Boltzmann distribution, but an arbitrary isotropic function. 

Specification of a set of macroscopic variables obeying an autonomous evolution in agreement 
with experiment is the aim of any non-equilibrium description of a system. The example we consider 
here is single-phase compressible fluid mechanics, where the average velocity, density and energy 
of a vast (~ 10 23 ) number of atoms or molecules moving in D dimensions obey a closed set of 
equations for D + 2 macroscopic variables. 

Boltzmann provided us with the first concrete steps towards a general method for finding such 
closed descriptions. The kinetic equation which bears his name is an approximation to the first 
equation in the BBGKY hierarchy. In order to obtain a closed equation for the single particle 
distribution function, Boltzmann made his famous stosszahlansatz, or molecular chaos approxima- 
tion, replacing the two-particle distribution function by a product of one-particle distributions. 
We may regard the stosszahlansatz not as a restriction on the behaviour of microscopic dynamics, 
with which it is in direct contradiction, but as a restriction on the types of macroscopic variables 
for which we are seeking autonomous descriptions. That is, we seek autonomous descriptions for 
macroscopic variables depending on the single particle distribution function. 
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Having obtained the Boltzmann equation, the Chapman Enskog expansion gives us a method 
for finding the moment equations [3JG3. The moment equations, as we shall see, are obtained as 
a hierarchy (in this case referred to as the super-Burnett hierarchy). The lowest order equations 
presume that an exact local equilibrium exists everywhere, and the resulting equations are the 
inviscid Euler equations. At the next level one introduces small deviations from equilibrium, and 
the resulting moment equations are the Navier Stokes equations. At each successive level in this 
hierarchy one expresses higher order corrections to the local equilibrium distribution in terms of 
gradients of that local equilibrium. 

As emphasized by Chapman and Cowling [I] , the division of the distribution into equilibrium 
and non-equilibrium parts is non-unique. It is therefore unclear what range of validity a given 
moment equation in the super-Burnett hierarchy has. There are known situations of consider- 
able practical interest in which the Navier Stokes equations must be replaced by a super-Burnett 
scheme In general it is necessary to validate any given scheme against some more fundamental 
approach such as direct simulation Monte Carlo 0. Furthermore the entire edifice crumbles if 
we do not know what equilibrium distribution the fluid will tend to. This has been a particular 
stumbling block when one wishes to extend the Chapman Enskog treatment to interacting fluids 
where one does not in general know the equilibrium distribution. 

However, even in the case of non-interacting fluids an interesting case presents itself. Tsallis 
and others have proposed various generalizations of Boltzmann-Gibbs thermostatistics, and experi- 
mental studies have verified that such statistics arise in nature |Sj ■ Recent work by Beck and Cohen 
has placed all such generalized statistics schemes on a common footing by utilising the physically 
appealing notion of fluctuating intensive parameters j^j. 

From an equilibrium point of view, such thermostatistics are somewhat controversial. It is 
unclear under what circumstances a system's equilibrium will deviate from the Boltzmann-Gibbs 
form. Although there are many systems known to so deviate, without a theoretical connection 
between a systems microscopic properties and its deviations from Boltzmann-Gibbs statistics, such 
distributions remain unsatisfactorily phenomenological. Away from equilibrium, the picture is quite 
different. Boltzmann's treatment of irreversible processes does not specify a single functional form 
for the entropy out of equilibrium. His .ff-function is equal to the Gibbs entropy only for an ideal 
gas in equilibrium. In any other circumstances these quantities are not equal JOJ- In general, 
the entropy away from equilibrium is only defined with respect to a set of macroscopic variables 
undergoing an autonomous evolution. Consideration of the widest possible range of distributions 
and entropies seems only prudent in this context. 

We have emphasized that it is the existence of autonomous evolution equations for macro- 
scopic parameters which defines the subject of non-equilibrium thermodynamics. Single phase 
fluid mechanics is the canonical example of such an autonomous description, and therefore how 
this description depends on the details of the local distribution is of some interest. In particular, 
under what circumstances does the hydrodynamic description cease to be autonomous? Conversely, 
determining the dependence of the hydrodynamic equations on the details of the distribution may 
provide experimentalists with the capability to determine the thermostatistics of a system by purely 
hydrodynamic experiments. 

The Boltzmann equation, H-theorem and Chapman-Enskog expansion have previously been 
studied in the context of Tsallis thermostatistics, and a set of Navier-Stokes equations have been 
derived for the Tsallis generalized statistics |2l llll l5]. Boghosian found that the incompressible limit 
of the Navier-Stokes equations is independent of the parameter q characterizing the Tsallis distri- 
bution, whereas the energy equation in the compressible case contains an additional g-dependent 
contribution to the heat flux This result was confirmed in and the additional term is the 
only further term allowed by the Curie principle |12j . 
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We begin in Section |2 by constructing the well known Chapman-Enskog expansion for the 
Boltzmann equation with a BGK collision operator. We then write the first and second order 
solutions in terms of moments of the distribution function. In Section |3] we introduce the "flux- 
field" expansion defined in ^3] an d the recursion relation proved in PP, and use this recursion to 
evaluate the moments, thereby obtaining the hydrodynamic equations for an arbitrary isotropic 
local equilibrium, in terms of the low order moments of that distribution. We close the paper with 
some further discussion and directions for future work. Throughout the following we choose units 
such that the sound speed c = 1 and the mass of the molecules composing the fluid m = 1. In 
particular this means that the hydrodynamic velocity u is given in units of the Mach number. 



2 Chapman - Enskog Expansion 

We begin with the BGK equation in D dimensions: 



d(f) 



■ ±(*-*W>). 



(1) 



Where cf> is the single particle distribution function, ^>(°) is an arbitrary isotropic function, v is the 
molecular velocity, r is a parameter (the collision time) and t is time. We define the operator 
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in terms of which we may write the formal solution of Q: 

(f) = (i+£>)-y o) 



(2) 



(3) 



We now introduce a hierarchy of separated length and time scales, such that we write the operator 
V: 



n=l n=l \ n / 



(4) 



Usually in the Chapman Enskog expansion one considers that relaxation to a local equilibrium 
happens on timescale t\, and slower kinetic and hydrodynamic effects occur on progressively longer 
timescales. We wish to consider a slightly more general situation where the distribution relaxes to 
a locally isotropic form on timescale t\. Subsituting the above in our formal solution (j2J) we obtain: 



i + £(_i)p r p(£; £ »z> n 

P =i 



\n=l 



Expanding our distribution about the isotropic part: 



+ + ■ ■ ■ 

yields a hierarchy of equations for the successive deviations from equilibrium cff 1 : 

0(2) = - T (V 2 - tVD^ 

0(3) = -t(V 3 - 2t£>i£> 2 + T 2 Vf)(p( 0) 



(5) 



(6) 



(7) 
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If we additionally specify that the deviations from equilibrium are not to affect the equilibrium 
values of the moments of the distribution we may obtain a hierarchy of conservation laws for those 
moments. At order e we obtain the continuity and Eulers equations for the fluid, at order e 2 
we obtain the Navier-Stokes equations for viscous compressible flow and at subsequent orders we 
obtain the Burnett equations and super-Burnett equations. In the present paper we consider only 
the first two orders in this hierarchy, but note in passing that the techniques employed here are 
equally applicable to obtaining such higher order hydrodynamic schemes. 



2.1 First order solution 

Multiplying the first order equation by the mass, momentum and kinetic energy of the fluid particles 
and integrating over velocity gives: 



^) d D v-rV- J vv ) <j&)d D v (8) 

2 / \ 2 / \~ 

We define the mass density p, momentum density pu, internal energy p£: 



p= J (fPd?V 

pu= v(j)°d 3 v (9) 



p£ = \J v 2 <f>°d\ 
and the inviscid pressure tensor and heat flux pj^: 

p j(0) = J ^0(O) vd 3 

and hence obtain the continuity and Eulers equations for the fluid. 

^ + V.(pu) = 



dti 



d ± + V>J<°>) = 



2.2 Second order solution 



(10) 



+ V-P u = (11) 



We first simplify the equation for the second order deviation from the equilibrium using the first 
order equation: 

^ = -r(V 2 - rVl)^ 
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Multiplying by the moments and integrating over velocities yields: 




0(2) v \d D v = -t4~ / v ^ 0) d D v - tV 




d 

r di 



(13) 



If we define the first order corrections to the pressure tensor and heat flux due to deviations from 
equilibrium: 



p3 W = J ^(i) vrf 



(14) 



We obtain the viscous Navier-Stokes equations: 



f + V.(pu)=0 



2 



d P u + V .( V (°) +V W) =0 (15) 



dt 2 

*g- + v • + ^JW) = 

The expressions for the corrections to the pressure tensor and heat flux may be expressed in terms 
of gradients of the local equilibria by using the first order solution to our Boltzmann equation. Of 
course, the above equations will only be a useful fluid dynamical scheme if the pressure tensor and 
heat flux may be expressed entirely in terms of other macroscopic quantities. It is to this issue that 
we turn in the remainder of the paper. 



3 Moments of the distribution function 

In this section we evaluate all moments of the distribution function defined above. The Chapman- 
Enskog expansion enables us to express deviations of the distribution function from equilibrium in 
terms of gradients of the local equilibrium distribution. All moments may therefore be evaluated 
in terms of the function cjy^' . As we shall see, the conditions of isotropy and Gallilean invariance, 
which this function must satisfy, are sufficient to enable us to evaluate all moments of this function 
in terms of the lowest order moments. These lowest order moments are precisely the hydrodynamic 
variables, and so we obtain an autonomous macroscopic description under very mild assumptions 
about the form of the local equilibrium distribution. 

Our results will follow from a remarkable property of the Taylor expansion of isotropic functions, 
which is proved in pQ. We first introduce some notation. We denote the mth rank outer product 
of a vector v with itself by: 

m 

V = WW • • • vw 

* , ' 

771 times 




5 



similarly we denote the rath rank outer product of gradient operators by: 

v = vvv--- vv 

v . ' 

m times 

this operator is equal to the familiar Hessian operator for m = 2 and is the higher rank generalization 
of the Hessian for ra > 2. 

We combine this notation for the outer product of vectors with a similar notation for the inner 
product (or more generally a p-fold contraction operation): 



r—p q—p 



aQ(g)b = (a bf (g)a(g)b 



(16) 



The above signifies the p-fold contraction of a qth rank tensor (in fact a g-adic) on a rth rank tensor 
(an r-adic). 

We now write the usual Taylor expansion of a function of a vector: 



/(x + a) = 2>-V)V(x) 

n=0 

We write this in the notation defined above as: 

oo n n i- n 

/(x + a) = £(g)aQ (g)V/(x) 

n=0 



(17) 



(18) 



While it may appear that the purpose of this notation is to take the simple, well known expres- 
sion for the Taylor expansion and convert it to a complicated and unwieldy one, we are not being 
deliberately obscurantist. The point of this notation is to associate with each order of the Taylor 
expansion the tensor: 



®V/(x) 



(19) 



We are interested in isotropic functions of vectors, that is, functions only of the modulus of the 
vector. As we shall see in a moment, the tensors T n associated with the Taylor expansion of such 
isotropic functions take a particularly simple form. 

We follow the notation of and define the nth-rank completely symmetric tensor kernel 



*n(v)= £ 



m=\n/2\ 



4>m(v) 

(n — m)\ 



per 



'2m— n 



(20) 



where "per" indicates a summation over all distinct permutations of indices, and where we have 
defined the following functions related to the derivatives of 4>(v), 



(v 



I d 

v dv 



[v . 



(21) 



and state the following property: 



V/C n (v) = /C„ +m (v) 



(22) 
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The first few such kernels are given in appendix 

We may write the Mach number expansion of our distribution as a sum of complete contractions 
of totally symmetric nth rank tensors: 

oo _. ■- n -inn 

0(|v-u|) = ^(-l) n - (g)V0(|v|) Q(g)u (23) 

n=0 ' L 

Using our tensor kernels defined above and the identity (|22[) we identify /Co(v) = 0(|v|) and write 

0(|v - u|) = E(-l)"^(v) O ® u (24) 

n=0 

For an isotropic function the tensors T n appearing in the Taylor expansion are exactly the com- 
pletely isotropic tensors given by (|2Uj) . If one compares this expression for the Taylor expansion of 
a function of the modulus of a vector with that given in |13j the economy of this notation becomes 
apparent. 

We now define the moments of the distribution via the two sets of integrals: 

„ m 

C = j v K n (v)<g)vd D v. (25) 

and 

1 r m 

J n=~ 2 \ 1C n (v)v 2 (g)vd D v (26) 

These integrals have the property I™ = if n > m and J™ = 0ifn>m + 2. Furthermore both 
integrals vanish if m + n is odd. This implies that for each moment of the distribution function 
only a finite number of terms in the Mach number expansion is necessary to exactly compute 
that moment. Furthermore, as a consequence of l)22D the integrals above satisfy a relatively simple 
recursion relation which means that they may all be computed in terms of the lowest order integrals. 

We first consider the integrals I™. First note that the integrand is odd for m + n odd, and so I™ 
vanishes for m + n odd. The asymptotic behaviour of our tensors is governed by the set of functions 
4> n . The kernel /C n (v) behaves asymptotically like v n (j) n . We require /Co(v) to be normalisable, that 
is the integral: 

r K^)d D v, (27) 



must converge. This means that if /Cq(v) is algebraic it must behave asymptotically like ~ v~ p 
where p > D — 1. The functions (f> n behave asymptotically like u - (? 5 + 2n ) ) and hence the kernel 
/C n (v) behaves asymptotically like v~ q where q = (p + n) and q > n + D — 1. This implies that the 
integrals: 

J K n (v)(g)vd D v, (28) 

converge if m < q, i.e. ifm < n + D — 1. The first two integrals are determined by the fact that 
/Co is the isotropic part of the comoving distribution function: 

i (29) 
'o=0 
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We define the temperature as proportional to the average of the particles kinetic energy, where the 
proportionality is chosen in accordance with the equipartition theorem. 



pDkT ^ 1 



< l mv 2 >=\J K (v)v D+1 dv Jd D n = ^J K (v)v D+1 dv 



(30) 



Where Sd is the surface area of the unit sphere in D dimensions. This definition of temperature 
is microscopic, and hence valid out of equilibrium. This will correspond with the macroscopic 
equilibrium thermodynamic temperature if and only if the comoving distribution is Maxwellian. 
We impose no such restriction. Introducing the unit vector e„ in the v direction we write our first 
nontrivial tensor moment in terms of the temperature: 



J2 = J Kv{v)v D+l dv j e v e v d D n = ^1 j K^)v D+1 dv 



Where we have used: 



Giving: 



j e v e v d D U = ^1 



il = pkTl 



We can now proceed to the general case. We state the following theorem: 



V (x) v = per 



(1) 



vper 



fm-2 



(31) 

(32) 
(33) 

(34) 



where by definition any outer tensor product of negative rank is zero. The proof of this theorem is 
given in (J. We utilise this below: 



V 



£ n _i(v)(S?)v 



/C n _i(v)(S?)v 



: n _!(v) }d D v 

m— 1 



dV 



J K, n -i{y) | per 



vper 



fm-1 



v (8)1 



(35) 



d D v 



The boundary term vanishes for m < D + n — 1 and we obtain the slightly awkward recursion 
relation: 



j /C n _i(v) 



per 



'm—l 



vper 



'm-2 



d D v 



(36) 



We write the recursion relation in terms of tensors per[m, 1], which are defined and discussed in 
detail in appendix 



C = - / X^_i(v)pCT[m - 1, l]d D v 



(37) 



We immediately obtain two results of importance. Firstly we note that I® = for n > 0. This 
property guarantees that the higher orders in the Mach number expansion of our equilibrium 
distribution do not contribute to the value of the density. Secondly, using the recursion relation it 
is clear that I® = for n > implies the general result that I™ = for n > m. 
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We now consider the integrals J™, which, by an extension of the argument given above for /, 
are convergent for m < n + D — 3. Write: 



m 
n i 



/C n _i(v)v 2 



mm m 

/C„-i(v)V( W 2 ) (9) v + K n {v)v 2 (9) v + /C n _i(v) W 2 V (9) v (38) 



giving: 



1 



^ 2 /c n _!(v; 



dV 



71-1 



771+1 

(v) (9) vd D W 



(^ 2 )/C n _ 1 (v)V 



Introducing our tensors per[m, 1] we obtain: 

<C = -C-l 1 ~\j (v 2 )IC n ^(v)p^[m - 1, l]d D V 



(39) 



(40) 



The temperature is defined via Jq = DpkT/2. The recursion above implies J° = for n > 2, and 
hence J™ = for n > m + 2. Evaluation of the first integrals gives: 



J° = 




= P i 






-ug 


[«^3 = 







pkT(D + 2) 



2 2 



(41) 



Where fiyjy is the completely isotropic fourth rank tensor defined in appendix [Bj We also require 
the integrals Jq , Jf , J| . The integral Jq does not follow from any of the integrals considered so 
far. This is because it is proportional to the fourth moment of the distribution function. For 
the Maxwell-Boltzmann distribution the fourth moment may be expressed in terms of the second 
moments of the distribution. For the quite general distributions considered here this may no longer 
be true, and we introduce a new scalar variable R: 



R 



1 



v A ^)d D v 



D+3 6dv / dU 



D(D + 2) J " v- ~ D{D + 2) 
ci 

We may now write Jq in terms of this new scalar variable: 



Sd 



D{D + 2) 



,-D+3 



4>dv (42) 



Where we have chosen the definition so that the related integral In is independent of dimension. 



7 2 



v 2 w(pd D v 



— T^Sn 
2D 



,D+3 



(frdv 



v 4 (f>d D v 
(D + 2 



v D+3 cf)dv J e v e v 



(43) 



Rl 



All integrals may now be computed in terms of p, T, R and u. The remaining higher order 
integrals are used below and computed in detail in appendix[Bj We begin with the inviscid pressure 
tensor V°: 



pi 



vv 



<P°o)d 3 



uu 



pkTl + puu 



(44) 
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the energy density and heat flux: 

H 2 
= P 
2 
1 

2 



v 2 <j>°d 3 v = J° + ^ : Jl 



pi 



DkT + u 2 ) 



j 2 ^°)d 3 



t1 UUU o t1 

-u • J{ - -5J- d J3 1 



(45) 



3! 



pA;T(L> + 2) pu 2 

u H u 

2 2 



Note that the internal energy is defined as: 



pi 



pDkT 



(46) 



and hence the ideal gas equation of state for the scalar pressure P is independent of the form of 
the distribution function when written in terms of the internal energy: 



P =D pl 

The first correction is to the pressure tensor: 

Substituting from above and using the integrals computed in appendix iBl 



(47) 



(48) 



dh 



+ rV 



u • If + 



3 , UUU n 3 r3 



(49) 



V • (pkTu) + V {pkTu) + [V (pkTu)] T + V • (puuu) 



The first order correction to the heat flux is then: 



pj 



(i) 



d_ 

dh 
d_ 

dh 



u-Jl + 



UUU 

~3! 
tV • Q 



Q 3 4 



rV 



2 UU 2 UUUU 2 

Jo + '■ J 2 "I : J A 



(50) 



Where we have introduced the heat flux tensor for brevity. This quantity may be written out 
explicitly as: 



Q 



(51) 



The only step which remains is to express the time derivatives with respect to h in terms of spatial 
derivatives using the first order equations of motion. After some laborious but straightforward 
manipulations we obtain: 



77(1) 



M[V(u) + [V (u)] t ] + A1V-u} 



(52) 
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Where we have introduced the shear viscosity: 

H = t pkT (53) 

and the bulk viscosity 

A = -T^pkT (54) 
We also obtain the first order correction to the heat flux: 

pjW = -u- ^i[Vu + (Vu) r ] + AlV-u^ +t ^ D ^ 2 ^ kT\7(pkT) - V(pR)^ (55) 
If we now define: 

6 = R - {kTf (56) 

We may write: 

pjW = - u • U[Vu + (Vu) T ] + A1V • u ) - T ^ D + 2 H pkTV{kT) + V{p6)\ 

V . ( 57 ) 

-u • f /i[Vu + (Vu) T ] + A1V • u ) - kVz - RV(p9) 

Where we have defined the thermal conductivity: 

D + 2 

k = -^—rpi (58) 

and the anomalous transport coefficient: 

D + 2 

R = -^-r (59) 

These transport coefficients have been defined such that the anomalous transport term vanishes 
when the quantity 9, vanishes. This is the case if 4>^ is a local Maxwellian, and if 4>^ is a Tsallis 
distribution these results are identical to those obtained in [2]. 



4 Conclusions 

We have shown that the Navier-Stokes equations for incompressible flow are independent of the 
detailed form of the local distribution function, provided only that the distribution may be written 
as an isotropic function of velocity with small anisotropic corrections. The energy equation in a 
suitably generalized form is also independent, acquiring only one extra term for non-Maxwellian 
distributions. This extra term is the only additional term allowed by the Curie principle. 

The original motivation of Boghosians derivation of hydrodynamic equations for systems with 
a local Tsallis equilibria was the potential for determining a systems thermostatistics by perform- 
ing entirely hydrodynamic experiments. The present work extends this hope to the full range of 
non-Maxwellian distributions currently being considered, with the caveat that two distributions 
whose fourth velocity moments have the same dependendence on their lower order moments will 
be empirically equivalent in this regard. 

We expect that the techniques employed in this paper will be of further assistance in the deriva- 
tion and elucidation of more complex hydrodynamic schemes by the Chapman-Enskog method. In 
particular, obtaining such schemes for interacting fluids is fraught with difficulty, and usually re- 
quires one to make some ansatz for the equilibrium distribution of the fluid |14| I15| I17j . One 
might hope that by using methods which do not require specification of the detailed form of the 
local equilibrium, one might place such interacting fluid models on a firmer theoretical foundation. 
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A Completely symmetric tensor kernels 

The first few tensor kernels as denned in eq. (|2Uj) are: 

[£i(y)li = Mv)yi 
[JQi(y)]ij = <f>i(v)8ij + <h(v)viVj 
[^3(y)]jj fc = 4>2(y)(yi s jk + yj5ik + yAj) + 4>3(y)yiyjyk 
[^(y)] ijk i = <fa(y) (Gijhi + $ik6ji + SuSjk) /2 + 

<fe(y) {yiVjhi + yiykSji + y^jfc + 2/j2/fc<5ii + yjyiSik + ykV&j) + 
<f>4(y)yiVjykyi- 



B Moments of flux field expansions 



We first give the value of the following angular integral, which will be useful for the higher order 
integrals to follow: 



eiejekeid D Q 



D(D + 2 

Sd 
D(D + 2) 



■{Sijbki + 5ik5ji + SuSjk) 



(60) 



n 



ijkl 



Where we have defined the completely isotropic fourth rank tensor fiy-jy. We consider the integrals 
if, 1^, if , if. The integrals I®, if vanish by symmetry: 

l\ = - J /C (v)per[0, l]d D v = -pi 

and 

If = - J /C (v)per[2,l]d D « 

= — j 4>v D+1 dv(^j ekeid D Vt5ij + j ejeid D Q5ik + J ejekd D 0,Sij\ 
= -pkTft ijk i 

We now consider the integrals 7^, > ^2 > ^2 j ^2 > ^2 • The integrals 7^, 7f vanish by symmetry, and 7^ = 
for n > 0, 



hijk 



J Ki(y)per[l,l]d D i 



4>i{v)v h {vj5 ik + v k 5ij)d D v 



1 

v dv 



— j-v 2 v D l dv[5ik j ehejd D Q + 5ij / ehe^dP^l 



Sd I 4>v D 1 dv[ 6 ik 6 hj + dijdhk 



P\ 5ik5hj + 5ij5hk 
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and 



efhijk 



£i(v)per[3, l]d D v 

(pi(v)v e (viVjV k Sfa + ViV k V[5 jh + ViViVj5 hk + VjV k vi5 ih )d D v 
= - J ~~j~ v4vD ~ ldv (^ hl J eeeieje k d D Q + 5 jh J e e eie k eid D Vt 
+ $hk J e e eieiejd D Q + 5ih J e e eje k eid D Q 



We know the result of the integration over solid angle and integrate over the speed by parts: 

Sd 



efhijk 



D 



J (f)V D+1 dv yShl^eijk + ^jh^eikl + ^hk^eilj + $ih®>ejkl\ 



Substituting again from our definition of temperature we obtain: 

= pkTl 5hl^l e ijk + Sjh^eikl + $hk^eilj + &ihQejkl 
efhijk V 

Finally, we consider Jf , the final integral required in the analysis above. 



J 3 



/C 2 (v)per[2, \}d D v 



[A 



3 J ijklr 



I ijlra 



&kr, 



~P [SjnSimSkl + Sj m 5i n 5kl + 5j n 5u5km + fijl5in<5km + ^jm^il^kn + ^jl^im^kn] 



We now turn our attention to the integrals Jf, j|, which may be defined in terms of previously 
computed integrals: 



4 = -Ii-\ f v 2 Kt(v)pSr[l,l]d D v 



pkT 



^ijkl 4"~ 



2 v L /ih J 1 2 

(D + 2) 



($i k 8ji + SuSjf.) 



and finally Jf : 



J 



2] 

4 J ijklmn 



v 2 )C 3 (v)pSr[l,l]d D v 

l3\ ijklmn ["^3 J ijkm ^ n \.'^^\ijkn^^ rn 



= P [Sjn^imSkl + Sj m 5i n 5kl + Sj n 5u5km + bjl&inbkm + Sj m 5u8kn + 8jl5i m 5kn] 
\S^ijkm$ln ^ijkn^lm] 
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C Permutation operators 



We have introduced the permutation operator per via: 

n 

per 




per [n, m\ 



where "per" indicates a summation over all distinct permutations of indices. The first few such 
tensors are thus given by 

per [0, 1] = 5ij 

per [0, 2] = 8ijS k i + 6 ik Sji + Su5j k 
per [1,0] = yi 

per [1, 1] = yi5 jk + yjS ik + y k 5ij 
per [2, 0] = yiyj 

per [2, 1] = yiyjSki + yiy k 5ji + yiyi5 jk + yjy k 5u + yjyAk + VkVAj 
per [3, 0] = yiyjyk 
per [4,0] = yiyjy k yi. 

We only need to consider permutations of the form per[n, 1] for the evaluation of our recursion 
relation. These permutations are: 

per [0, 1] = 5ij 

per [1, 1] = yi5 jk + yj5 ik + y k Sij 

per [2, 1] = y { (yj5 k i + y k 5 fl + y t 5 jk ) + y k yi5ij + yjyi5 ik + yjy k 5 u 

per [3, 1] = y h (yiyj5 kl + yiy k 8ji + yiyiS jk + yjy k Su + yjyAk + y k yAj) + 

{yiyjVkhi + ViVkVi^jh + ViViv^hk + yjVkyAh) 

per [4, 1] = yfiyhViVjSki + yhUiVkSji + VhUiy^jk + yhViVk^n + VhVjyAk + 

yhVkvAj + yiVjVk^hi + yiVkVi^jh + ViViv^hk + ViVkvAh) + 
ViUjUkyiShf + VhVjykyiSif + yhViVkyiSjf + yhynjjyihf + yhyiyjykhf. 

We have positioned the brackets above to illustrate the pattern which occurs among these tensors. 
Let P(n) be the number of terms in the permutation per[n, 1]. In the permutation per[n + 1, 1] 
there will be P(n) terms in which the new index is attached to one of the y's. per[n + 1, 1] is an 
n + 2th rank tensor and so there are n + 2 distinct ways of attaching the new index to the kronecker 
delta. Therefore the number of terms in per[n + 1, 1] is P(n + 1) = P(n) + n + 2. 
Our second permutation operation is defined in terms of the above operation: 

per[n, 1] = per[n, 1] — yper[n — 1,1]. 

Given the above decomposition of the per operation, it is straightforward to write down the first 
few such tensors: 

per [0,l] y = 5 tj 

Per [1, l]ijk = Vjdik + yAj 

Pe r t 2 > Mijki = VkVAj + VjVAk + yjyk^u 

per [3, i\hijki = (ytyjykbhi + yiykyrfjh + yiyiyjShk + VjUkvAh) 

per [4, l] fhijk = yiyjy k yi$hf + yhyjykyAf + VhViykyiSjf + yhViVjyihf + yhyiyjykSif. 
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Let P(n) be the number of terms in per[n, 1]. Using our previous recursion relation for P(n) it is 
straightforward to see that P(n) = n + 1. 
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